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The phenomenon of prethermalization and the subsequent steps of thermalization 
are analyzed in the framework of the chiral quark model. We solve the quantum 
equations of motion of the field theory derived from the 2PI effective action and 
■ study the time scales of equilibration. We find that already after a 0.6 fm/c long 

period of time some equilibrium features appear, even though the system is still 
' far from equilibrium. This might be an ingredient for understanding the success 

of ideal hydrodynamic description. 
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1. Introduction 



One of the crucial issues in heavy ion physics is the early thermalization of 
the excited quark and gluon fields: The extreme success of hydrodynamic 
. calculations started earlier than I fm/c require the theoretical establishment 

' of the early existence of an equation of state. 1 

The proper description of equilibration is a longstanding challenge in 
particle physics. Kinetic theories are popular tools for studying the series of 
incoherent collisions. By neglecting the off-shell processes the easiest way is 
to follow the clastic collisions only, but these deliver thermalization times by 
orders of magnitudes larger than desired. Equilibration may be enhanced 
by inelastic scatterings so that it completes before chemical freeze-out 2 . 
This way, thermalization times reduce to 2-3 fm/c, which is still longer 
than the early use of flow equations would require 3 . 

Recent field theoretical developments enable the nonequilibrium treat- 
ment of quantum fields by performing a systematic ladder resummation 4 . 
The available numerical techniques and resources enable us to solve the 
time evolution of 3+1 dimensional quantum fields in a self-consistent ap- 
proximation scheme. These methods are not limited to scalar fields any 
more, nowadays the 3+1 dimensional chiral quark model (with two quark 
flavors and 4 scalars) may also be routinely solved 5 . 
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Figure 1. Left: The fermionic occupation number n^\t;p) for three different mo- 
mentum modes as a function of time. The evolution is shown for two different initial 
conditions with same energy density. The dynamics becomes rather quickly insensitive 
to the initial distributions (A) and (B) displayed in the insets - much before the modes 
settle to their final values. Right: The ratio of pressure and energy density as a func- 
tion of time. One observes that an approximately time-independent equilibrium form 
of the equation of state builds up very early. The inset demonstrates for two different 
couplings that this prcthcrmalization time is independent of the interaction details. The 
insensitivity of bulk quantities to the late-time thermalization scale is in sharp contrast 
to mode-by-mode quantities. The long-time behavior is shown on a logarithmic scale for 
t > 30m- 1 . 

Thermalization time is not a single number, of course. Different degrees 
of thermalization (LTE, kinetic and chemical equilibration) can have very 
different time scales ranging from the rapid convergence of bulk quantities 
to the slow equilibration of energy spectra 6 . 

The model we consider is built up from two massless fermionic degrees of 
freedom and four scalars interacting according to the following Lagrangian: 

L = + | [(da) 2 + (Ott) 2 ] + [a + i l5 r a TT a } i> ~ V(a 2 + tt 2 )}, 

where r a denote the standard Pauli matrices. We consider a quartic scalar 
self-interaction V(a 2 + tt 2 ) = m 2 (a 2 + tt 2 ) /2 + A (a 2 + tt 2 ) 2 /{A\N 2 ). The 
employed couplings are taken to be g = 1 and A = 1. Though these 
couplings are smaller than heavy ion applications would require, we will be 
able to draw conclusions for the physical values since our main result will 
turn out to be coupling independent. The smaller-than- natural coupling 
leads to a spectacular separation of the various time scales. 

We use the two-particle irreducible (2PI) effective action to two-loop or- 
der, which includes direct scattering as well as off-shell and memory effects 5 . 
We numerically solve the noncquilibrium gap equation and obtain the evo- 
lution of the propagator G(ti,t2, \x\ — X2\) for each degree of freedom. To 
simplify our equations we assume homogeneity and isotropy in space. 
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Figure 2. Left: The solid line shows the kinetic temperature T^ n (t) as a function of 
time. In contrast to the mode-by-mode quantities, the kinetic temperature, which in- 
cludes all momentum modes, exhibits the phenomenon of prethcrmalization. Also shown 
is the ratio of the scalar and the fermion chemical temperatures, T f (t)/T / (t), for 

1 ' enem^ '< enem v 

different values of the coupling g and T eq . One observes hat chemical equilibration does 
not happen on the prethcrmalization time scale. Right: On the damping time scale a 
non-equilibrium variant of the fluctuation-dissipation relation is dynamically realized. 
This gives a particle number definition even for strongly coupled systems. The grey 
bands demonstrate the particle numbers and the consistency of their determination. 



2. Time scales of thermalization 

Thermalization is loss of initial information. Our dynamics is conserva- 
tive, still, virtually, we loose information if we do not read out all the 
internal variables but only the time-local propagator G(ti ~ t 2 ,\xi — x 2 \). 
This time-local propagator stores the information on the particle number 
distribution 5 . The evolution of these occupation numbers for three differ- 
ent momenta is shown in Fig. lleft. We also show the convergence for two 
different initial spectra with equal energy density. Although the details of 
the initial distribution are mainly washed out in the rather early evolution 
(15|m| _1 ), the spectrum is still far from equilibrium. It takes by an order of 
magnitude longer to approach the thermal spectrum (for small coupling) . 

We may split the observed time scales to three epoches. The separation 
of these scales we find for perturbative couplings only, in a strong coupling 
application the first two might coincide. (The time values correspond to 
the parameters used in the presented plots.) 

(1) Prethermalization — loss of coherence (t < 5|m| _1 ) 

Prethcrmalization is a universal far-from-equilibrium phe- 
nomenon, its time scale is given by the relevant scales of the dy- 
namics (temperature, mass) but it is independent of the coupling, 
as shown in Fig. 1 right. 
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Contrary to the mode-by-mode occupation numbers (Fig. lleft) 
the bulk, averaged quantities are often insensitive to the details of 
the spectrum as it is also observed in classical field theory 7 . 

We find that this short time is enough for very rapid establish- 
ment of an almost constant ratio of pressure over energy density 
(Fig. lright), as well as a kinetic temperature based on average 
kinetic energy (Fig. 21eft, for definitions see Ref. 6 ). These are im- 
portant ingredients of the early equation of state. This level of 
equilibration is one of the important preconditions of the use of 
hydrodynamic equations. 

(2) Damping — decay of individual modes (t < 15|m| _1 ) 

The individual excitations are washed out so that a smooth but 
still far-from-continuum spectrum is formed. We find that this re- 
laxation time agrees with the inverse width (7(2?; t)) of the nonequi- 
librium spectral function 5 . 

We know that the damping time (obtained from the imaginary 
part of the self energy or from the width of the spectral function) 
is the decay rate of an individual signal excited over the ensemble 8 . 
This relaxation time may be much shorter than thcrmalization 
time 9 . In the leading order perturbation theory relaxation time 
is proportional to the coupling-squared. We also see that this re- 
laxation time scale is not a subject of spectacular evolution. 

Within this time scale the generically independent symmetric 
propagator (F) and spectral function (p) evolve to fulfill a relation 
(in Wigner coordinates): F(u),p; t) — ±(n(u>; t) + ^)p(w,p; t), which 
is known as fluctuation -dissipation relation in equilibrium. Here 
the nontrivial issue is the momentum (p) independence of n(u);t). 
This enables us to translate it to n(p; t) using the dispersion rela- 
tion defined by the peak of the spectral function. The existence 
of the relation above provides a particle number definition n(p;t) 
not tainted by quasi-particle picture. 10 Fig 2right shows by shaded 
bands how consistently the particle number n(p;t) is determined 
from the evolution of different k modes with Uk ~ w p ± 7 P . 

(3) Thermalization — loss of spectral information (t < 500|m| _1 ) 

Kinetic equilibration is a result of the collective evolution of the 
particle spectra. At the same time the particle abundances in the 
various species arc equilibrated as well. Whether this longer epoch 
fits into the short time range before chemical freeze-out is still a 
matter of discussion. 
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3. Discussion 

Direct heavy ion applications require higher couplings than in the example 
above. Tuning the coupling to g = 2.5 — 3 one observes that both relaxation 
and thermalization time shrinks considerably, approaching the prethermal- 
ization time. The elementary oscillators are then over-damped, and the 
estimated prcthcrmalization time will also give a hint for the damping time 
scale. With these couplings one can easily obtain thermal spectra in the 
fcrmionic sector within 1 — 3 fm/c. 

We find that prcthermalization time can be estimated by the inverse av- 
erage kinetic temperature 6 . In particular, we obtained the following quan- 
titative relation: Tt pt ~ 2 — 2.5. If one attempts to use the concept of 
prethermalization for the case of heavy ion collisions one has to replace 
this temperature by the relevant scale of the dynamics. This we assume 
to be the saturation scale: Q s « 1 GcV at RHIC energies. Estimating the 
prethermalization time by 3 times 1/Q S one has: i pr cthermaiization ~ 0.6fm/c 
This is the very close to the instant that is usually taken as initial time for 
the hydrodynamic evolution 11,12 . 

If the complete local equilibration cannot be established, it is worth 
thinking if this partial equilibration can already establish the use of flow 
equations. We do believe that the study of partial equilibration can help us 
overcome the discrepancy between the long thermalization time predicted 
by kinetic theories and the short times used by hydrodynamic models. 

The author acknowledges the fruitful discussions and collaborations on 
this work with J. Berges, J. Serreau and C. Wetterich. 
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